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ABSTRACT 

Over the past thirty years, laboratory experiments have been devised to develop fatigue crack growth rate data that is 
representative of the material response. The crack growth rate data generated in the laboratory is then used to predict the 
safe operating envelope of a structure. The ability to interrelate laboratory data and structural response is called similitude. In 
essence, a nondimensional term, called the stress intensity factor, was developed that includes the applied stresses, crack 
size and geometric configuration. The stress intensity factor is then directly related to the rate at which cracks propagate in a 
material, resulting in the material property of fatigue crack growth response. Standardized specimen configurations and 
experimental procedures have been developed for laboratory testing to generate crack growth rate data that supports 
similitude of the stress intensity factor solution. In this paper, the authors present laboratory fatigue crack growth rate test data 
and finite element analyses that show similitude between standard specimen configurations tested using the constant stress 
ratio test method is unobtainable. 


Introduction 

One of the responses of a material to extreme forces, such as stress, temperature, etc., is to crack. A crack appears when the 
material reaches a limit in its capability to absorb damage and fails. Sometimes, a crack will grow under a periodically applied 
condition, such as cyclic loading, that are well below the stresses required to fail the material, denoted fatigue crack growth. 
Over the past thirty years, laboratory experiments have been devised to develop fatigue crack growth rate data that is 
representative of the material response. The crack growth rate data generated in the laboratory is then used to predict the 
safe operating envelope of a structure. The ability to interrelate laboratory data and structural response is called similitude. In 
essence, a nondimensional term, called the stress intensity factor, was developed that includes the applied stresses, crack 
size and geometric configuration. The stress intensity factor is then directly related to the rate at which cracks propagate in a 
material, resulting in the material property of fatigue crack growth response. 

Standardized specimen configurations have been developed for laboratory testing to generate crack growth rate data that 
supports similitude of the stress intensity factor solution. Recent research into fatigue crack growth rate data has exposed 
some limitations in the testing standards [1, 2]. The typical approach to generate crack growth rate data at a specific stress- 
ratio is accomplished by reducing both the maximum and minimum applied load at a specific rate until the crack arrests. 
Several researchers [3 - 6] have all postulated that a constant stress ratio load reduction method produces results that are not 
representative of the material behavior but an artifact of the laboratory environment. Furthermore, Liknes and Stephens [7] 
and Garr and Hresko [8] have suggested that specimen configurations contained within the standards can have an effect on 
threshold by influencing extrinsic closure affects. The middle crack tension, compact tension and eccentrically loaded edge- 
crack tension specimen configurations are widely used for generating fatigue crack growth rate data [9]. In this paper, the 
authors present laboratory fatigue crack growth rate data generated using the standard constant stress ratio load reduction 
method that is specimen dependent and does not support similitude. 

Experimental Method 

Fatigue crack growth rate data for D6AC steel is presented in this paper that is dependent on specimen configuration. The 
steel alloy D6AC was chosen for study because it is used extensively by the aerospace industry [9]. Constant stress ratio, R, 
of 0.1 and maximum stress intensity factor range, K ma x, of 22 MPa m 1 ^ fatigue crack growth rate data was generated using 
four different compact tension, C(T), specimen configurations (1/1/ = 76.2 mm, B = 12.7 mm, W = 76.2 mm, B = 5.08 mm, W = 



50.8 mm, B = 5.08 mm, and W = 25.4 mm, B = 5.08 mm), one eccentrically-loaded side edge tension, ESE(T), specimen 
configuration (W = 38.1 mm, B = 5.08 mm), and two middle tension, M(T), specimen configurations (1/1/= 76.2 mm, B= 12.7 
mm, and 1/1/ = 76.2 mm, B = 5.08 mm). The laboratory environment, denoted “room temp, lab air” varies between room 
temperatures of 20° to 30° C, with a mean of 23° C, and the laboratory air relative humidity varies between 24 to 38%, with a 
mean of 28%. All experiments were conducted in loading frames of 89 to 222 kilonewton capacity. Load range cards were 
used to improve resolution for lower applied forces to within 2% of the expected force. The tests were conducted in force 
control at a frequency of 20 Hz using an analog controller. A computer control system [10] was used for K control of the C(T) 
specimens to adjust the forces based on compliance readings from strain or displacement gages and user input. Visual 
measurements of the crack length were taken throughout the test to verify the compliance based crack length used to control 
the test system. The M(T) specimens were installed in hydraulically actuated wedge grips and force was shed manually using 
the step method described in ASTM E 647. The C(T) and ESE(T) specimens were installed in sized clevises for each width 
tested and force was shed using K control at a rate of -80/m. A minimum of one repeat was performed for all test conditions 
presented, on a separate specimen. The scatter was less than 10% on AK for any given growth rate, da/dA/for all conditions 
tested. 


Computational Method 

Three-dimensional finite element analyses of both the C(T) and M(T) were analyzed using the ZIP3D [11] finite element code. 
ZIP3D is an elastic-plastic material, non-linear finite element software with capabilities to carryout fatigue and fracture analysis. 
A typical detail of a ZIP3D finite element model along with a schematic of the C(T) specimen are shown in Figure 1 . The C(T) 
was modeled with a thickness B = 12.7 mm and width W = 50.8, and an initial crack length to width ratio, ai/W = 0.25. . The 
analysis has not been conducted on the tesed ductile steel, D6AC, but the trends in crack opening displacement presented 
herein should be similar. The material considered throughout the analyses was aluminum alloy 7075-T73. The material was 
assumed to have a bi-linear character with modulus E = 71 .0 GPa and flow stress o 0 = 400 MPa The von Mises yield criterion 
and the associated flow rule were used. Small deformation theory was employed. The load cycles were applied to simulate 
load reduction from as high a AK of 30 MPaVm all the way to 1 MPaVm under both constant R and constant K max load 
reduction procedures. The crack front and load were adjusted during each analysis step such that da = 10 microns (one 
element size). Complete information on these analyses can be found in Seshadri and Forth [12]. 
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Figure 1 — A typical ZIP3d finite element model of C(T) specimen. 

Constant Stress Ratio Load Reduction 

Crack growth rate data was generated using M(T) and C(T) specimens of varying thicknesses and widths at a stress ratio, R, 
of 0.1 . Figure 1 plots crack growth rate, da/dN, versus stress intensity factor range, AK, for a stress ratio, R = 0.1 , condition 
omitting the 50.8 and 25.4 mm wide C(T) specimen data. The differences in crack growth rate data between the specimen 
configurations presented in Figure 1 are large as the growth rate drops below 10" 8 meters/cycle. The divergence of the 
specimen data occurs because some of the specimens are developing environmental- and roughness-induced fretting-oxide 


debris crack closure that retards the crack [5]. Photographs of representative fracture surfaces for the 5.08 and 12.7 mm thick 
C(T) specimens are shown in Figure 3 to illustrate the development of fretting-oxide debris on the 12.7 mm thick C(T) 
specimen fracture surfaces (Figure 3a). The 5.08 mm thick C(T), as shown in Figure 3b, M(T) and ESE(T) specimens do not 
develop oxide debris on the fracture surface at the humidity levels seen in laboratory air [13]. A comparison of the local cyclic 
CTOD and crack closure levels for M(T) and C(T) specimens computed in the finite element analyses is shown in Figure 4. 
The black and red lines show variation in cyclic CTOD across the thickness with applied AK for M(T) and C(T) specimens 
respectively. The blue and green lines represent variation in closure level across the thickness for M(T) and C(T) specimens. 
From the beginning of the constant R load reduction procedure, the estimated cyclic CTOD values across the thickness for 
M(T) specimen are within the bounds of C(T) specimen. That means to say that at any applied AK during the constant R load 
reduction procedure, the estimated local cyclic CTOD at any location across the thickness of M(T) specimen is a few percent 
less than the corresponding C(T) specimen. Also, during the initial part of the constant R load reduction procedure, the M(T) 
specimen has higher local crack closure level across the thickness when compared to the C(T) specimen. However, once the 
applied AK reduces below 6 MPaVm, locally the crack remains fully open for both specimen types. In general, during the 
constant R load reduction procedure, both M(T) and C(T) specimens have similar local cyclic CTOD behavior. Therefore, the 
crack growth rate behavior shown in Figure 2 is not a function of plasticity-induced crack closure, but of a specimen 
configuration susceptibility to the development of fretting-oxide debris under the constant R= 0.1 load reduction test method. 
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Figure 2 — Crack growth rate versus stress intensity factor range at R = 0.1 for different specimen configurations. 



(a) 

Figure 3— Photographs of representative fracture surfaces for C(T) specimens in lab air: (a) 76.2 mm wide, 12.7 mm thick; (b) 

76.2 mm wide, 5.08 mm thick. 
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Figure 4— Comparison of variation in local ACTOD with applied constant-R load reduction procedure for C(T) and M(T) 

specimens. 


Constant K ma x Load Reduction 

Constant K max = 22 MPa m 1/2 tests were performed to develop crack growth rate data that has minimal influence from extrinsic 
effects, such as the environmental- and roughness-induced crack closure depicted in Figure 3a. The constant K max crack 
growth rate data plotted in Figure 5 shows that for any value of AK, there is little difference between any of specimen 
configurations. The lack of sensitivity of the K max test data to specimen configuration implies that the extrinsic effects, such as 
crack closure that can develop in a constant R = 0.1 load reduction test, have different effects on different specimen 
configurations. Comparison of local cyclic CTOD and crack closure levels for M(T) and C(T) specimens computed using the 
finite element method under the constant K max load reduction procedure are shown in Figure 6. Once again, the black and red 
lines show variation in cyclic CTOD across the thickness with applied AK for M(T) and C(T) specimens respectively. The 
closure levels are not shown because crack closure is approximately zero in a constant K max test. At any given applied AK 
during the load reduction procedure, the cyclic CTOD estimate for the C(T) is greater than the M(T) specimen. Plowever, the 
difference in the magnitude of the CTOD across the specimen thickness is small. Since the constant K max test procedure 
develops crack growth rate data that is independent of specimen configuration, it is representative of a material response, and 
provides a baseline of data to support similitude, i.e. the data can be used to predict the response of different structural 
configurations. 
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Figure 5 — Crack growth rate versus stress intensity factor range data for K ma x = 22 MPa m 1/2 in lab air. 
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Figure 6 — Comparison of variation in local ACTOD with applied constant-K ma x load reduction procedure for C(T) and M(T) 

specimens. 


Summary 

Fatigue crack growth rate data was generated and modeled using several different standard specimen configurations. The 
data was generated using two experimental procedures that have been standardized with the intent of generating fatigue crack 
growth rate data that supports the theory of similitude. In other words, the fatigue crack growth rate data generated with either 
the constant stress ratio or K max procedure is purported to be representative of a material response that can be applied to any 
structural analysis. However, tests conducted with the constant stress ratio load reduction procedure generated results that 
were dependent on the specimen configuration (Figure 2) contrary to the finite element model (Figure 4). The finite element 
model captures the effects of plasticity on crack growth and clearly shows in Figure 4 that crack growth rate is specimen 
independent. Therefore, the data presented in Figure 2 is principally affected by other forms of closure not modeled by the 
finite element analyses, such as environment and roughness. The data generated using the constant stress ratio load 
reduction test is only useful for understanding the behavior of a crack in a C(T), ESE(T) or M(T) specimen under the specific 
loading and environmental conditions. To extrapolate the fatigue crack growth response presented in Figure 2 to predict the 
response of a structural component would be unwise and potentially catastrophic. For it is unknown if the data presented in 
Figure 2 is conservative with respect to the response of a structure. In contrast, both the fatigue crack growth rate data (Figure 
5) and finite element analyses (Figure 6) generated using the constant K max test procedure were insensitive to the specimen 
configuration. Figure 7 is a plot of the fatigue crack growth threshold (where da/d/V ~ 10' 10 meters/cycle) versus stress ratio 
that clearly shows the threshold data generated using the constant R = 0.1 procedure varies from 3.5 to 9.76 MPa m 1/2 
dependent on specimen configuration, whereas the constant K max test procedure varies from 2.62 to 3.02 MPa m 1/2 . In 
summary, the constant stress ratio test data presented in this paper is not representative of similitude contrary to the finite 
element analyses and is therefore not useful in predicting the response of any other structural configuration. 
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